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Abstract. Consider the 3D Anderson model with a zero mean and bounded 

5— | ' i.i.d. random potential. Let A be the coupling constant measuring the strength 

t*H i of the disorder, and u(E) the self energy of the model at energy E. For any 

■^^ , e > and sufficiently small A, we derive almost sure localization in the band 

E < — ct(0) — A . In this energy region, we show that the typical correlation 

length ^e behaves roughly as 0((|E| — a(E))~ 1 ' 2 ), completing the argument, 

outlined in the unpublished work of T. Spencer I18| . 

PL,- 

<—| ! 1. Introduction, main result and steps of the proof 

+-> . 

1.1. Introduction. In this paper we want to carry out the program, sketched in 

the unpublished preprint of T. Spencer [T5], regarding the localization for the 3D 
l — '■ Anderson model in the so-called Lifshitz tails regime. 

The Anderson operator H* on the lattice 1? acts on the vector ip e Z 2 (Z 3 ) as: 

C& : (H^)(n) := -i(A^)(n) + W u {n)^{n) , (1.1) 

~f\ , where A denotes the discrete Laplace operator, 

(A^)(n) = Yl 1>(n + e) - 6ip(n) . 

eeZ3, |e|=l 

QQ \ We will assume throughout this paper 

Al The values of the random potential V u (-) are independent, identically dis- 
^- . tributed variables, with even, compactly supported, and bounded probabil- 

yj ' ity density p. 

r> A2 For any m £ N, E(V^ 2m (x)) < c with some constant c, and E (V%(x)) = 1. 

Let e(p) denote the dispersion law, associated with the Fourier transform of the 
Laplacian, (!FAf)(p) = —2e(p)f(p), where 



f(p) := (JT/)( P ) := ]T ^ P ' n f(n) , P G T 3 := [-1/2, 1/2] 
with its inverse 



n<El 



){n)= / cPpe a ^ n f(p). 

JT3 



One then computes 



3 

e(p) = 2^sin 2 (^-e Q ), (1.2) 

a=l 

where e a is a unit vector in the a direction. The spectrum of the unperturbed 
operator 7J° is absolutely continuous and consists of the interval [0,6]. 

l 
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In what follows we will denote by A(x, y) the kernel of the linear operator A 
acting on Z 2 (Z 3 ) (that is A(x,y) = (S y ,AS x ), where 5 X is an indicator function of 
the site x £ Z 3 , and (•, •) denotes the inner product of l 2 (Z 3 )). We will use the 
concise notation J in place of X- T 3-,fc whenever it is clear from the context that each 
of the k variables of integration is integrated out over a torus T 3 . 

We will investigate the properties of H* for a typical configuration oj in a weak 
disorder regime, namely at the energy range 

\ x2 f d 3 p 4 _ 

Aa,-A / -r^-A 
Jt3 e(p) 

with a = min{a; : x € support p} for any e > and A > being sufficiently 
smalo Most of the mathematical results on localization for operators with random 
potential in dimensions d > 1 have been derived using the multi-scale analysis 
introduced by Frohlich and Spencer [H] and by the fractional moment method of 
Aizenman and Molchanov [3J (we are going to use the latter approach). By now 
there exists extensive general literature on the Anderson localization problem, see 
for example [TH] and references therein. 

The quantity of the most interest is the typical asymptotic behavior of the so 
called Green function (also known as the two point correlation function, the prop- 
agator) 

R(x,y) = (HZ + E + ir,)- 1 ^^) 
in the limit rj \ 0. It plays a crucial role in determining, for instance, the conduc- 
tivity properties of the physical sample (whether it is an insulator or a conductor at 
a given energy band). On a mathematical level, investigation of the propagator can 
yield an insight on the typical spectrum of HQ, at the vicinity of —E. The Ander- 
son model (II. ip is characterized by the following dichotomy, [3] : Either the typical 
Green function R(x,y) decays at least exponentially fast when \x — y\ — ► oo, or it 
cannot decay faster than \x — y\~ 6 in a three dimensional case. The former behavior 
is called localization and necessitate that the spectrum is pure point almost surely 
in the vicinity of —E, [SJ. The localization region is naturally characterized by the 
so called correlation length (g, that is a typical length scale |a; — y\ at which R(x, y) 
starts to decay at least exponentially fast. We stress here the energy dependance 
of the correlation length, for it is going to play a role in our analysis. The con- 
sensus among the condensed matter physicists is that in 3D, in the weak disorder 
regime, there should be a spectral transition from point spectrum to continuous 
one. This phenomenon is known as the Anderson transition. The proof of such 
an actuality presents a great challenge in this subject. The region of the possible 
spectral transition is called the mobility edge. One can get a certain indication of 
the existence of the delocalized regime for the Anderson model by showing that 
there occurs a threshold energy Eq (which presumably coincides with the mobility 
edge of the problem) such that £e diverges as E approaches Eq. The present work 
can be seen as a step in this direction. 

The occurrence of localization at energies near the band edges at weak disorder 
is related to the rarefaction of low eigenvalues, and was already discussed in the 
physical literature by I. M. Lifshitz in 1964, see Section 3 in [13], and [H]. As far 
as the rigorous results are concerned, let us only mention the three closely related 
works: M. Aizenman [1] showed that the spectrum of H^ consists (almost surely) 



\a = inf(E : E G a(H^)) almost surely, see e.g. |19| for details. 
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of the localized eigenvalues at the energy range [Xa, Xa + X a ], with a — 5/4. This 
result was later improved by W-M. Wang [20] (a — 1), whose result was in turn 
enhanced by F. Klopp [TU], who extended the region of localization all the way up 
to the (negative) energies of order A 1+1 / 6 in 31?. In this work we push the upper 
bound of the localization region further up, to the value —CX 2 , and examine the 
behavior of the correlation length as a function of energy. 

1.2. Results. In order to formulate our main technical accomplishment we need 
to introduce some further notation. The self energy term <j(E), associated with 
H^, is given by the solution of the self-consistent equatioro 

a(E) = X 2 [ .. d l P -- . (1.3) 

J T 3 e(p) +E- a(E) y ' 

It is easy to check that <j(E) is positive and uniformly bounded by CX 2 with some 
constant C, provided that 

d 3 p 



E > En := X z , 

e(p) 

and E* := E — a(E) > for such values of E (the relevant properties of the solution 
of (|1.3[) are collected in Appendix [XJ Moreover, if 

E>E e (X):=X 2 f ^ + A 4 "% (1.4) 

J T3 e(p) 

then E* > CX 4 ~ e for an arbitrary small e and sufficiently small values of A. We 
therefore can define a (renormalized) free Green function 

R r (x, v ) = (-Ia + E- o{E) + zO)- 1 ^, y) (1.5) 

for every E in the above energy range. Let us also denote by 

R(x 1 y) = (H* + E + tO)- 1 (x,y) 1 (1.6) 

which is well defined a.s., [3]. 

The hallmark of localization is rapid decay of G(x, y) at energies in the spectrum 
of Hu, for the typical configuration uj. Rapid decay of the Green function is related 
to the non-spreading of wave packets supported in the corresponding energy regimes 
and various other manifestations of localization whose physical implications have 
been extensively studied in regards to the conductive properties of metals and in 
particular to the quantum Hall effect. Our main result (Theorems Q] and [5] below) 
establishes this behavior of the Green function at the band edges of the spectrum, 
by comparing it with the asymptotics of the free Green function R r (x,y). The 
behavior of the latter for \x — y\ ^> 1 is known [9]: 

Rr(x,y)- T ^— l e-^"^-y l. 
2x\x-y\ 

The implication is that the correlation length for the free Green function is (E*)~ x / 2 . 
Moreover, for the energy range (|1.4j) we have 

(| 2 ;- 2/ | + l)- 1 »A((| a ;-j/| + ir 1 /2 

whenever 

\x-y\<{E*)- 1 '\ 



Note that since we will be interested at negative energies, E in Hl,3t is assumed to be positive. 
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hence 

l*r(z,2/)| a » j, f-i-iwa ' ior\x-y\<(E*)- 1/2 , E>E e (X). 

With this estimate in mind, we present 

Theorem 1 (Local fractional moment bound). For H* as above, for any s < 1/2 
and e > 0, there exists Ao(e) such that for all X < Ao(e) and E > E e (X), one has a 
bound 

E \R(x,y) - Rr(x,y)\ s < 1 (s)- f, 1W2 . (1-7) 

{\x-y\ + l) 8 / 2 

with Ci(s) < oo, which holds for any pair {(x, y) E Z 3 x Z 3 : \x — y\ < (-B*) -1 ' 2 }. 

The estimate (jl.71) can be interpreted as follows: The typical correlation length 
S,e within this energy range cannot be smaller than 0((i?*) -1 / 2 ), in particular it 
grows as E* approaches zerqj The next proposition shows that £e cannot be of 
much greater scale either. 

Theorem 2 (Global fractional moment bound). For H^, as above, there exists 
Ao(e), so that for all A < Ao(e), s < 1/4, and E > E e (X), one has a bound 

E\R( Xl y)\ s < ^WeC^MB'JM (i. 8 ) 

A 

for all x, y e Z 3 . 

Let us list several known implications of the global fractional moment bound: 

i. Spectral localization (^\): The spectrum of H u within the interval (|1.4| is 
almost-surely of the pure-point type, and the corresponding eigenfunctions 
are exponentially localized. 

ii. Dynamical localization ([\\): Wave packets with energies in the specified 
range do not spread (and in particular the SULE condition of [TB] is met): 

\n V \{P {H ^ Et} e- UH ){x,y)\] < Ae'^^ , (1.9) 

\t£R J 

where Pu<a stands for the spectral projection of H on the energies below 
a. 
iii. Absence of level repulsion (|15j J. Minami has shown that (|1.8[) implies that 
in the range (|1.4[) the energy gaps have Poisson-type statistics. 

For energies E slightly above — Eq it is expected on the physical grounds (see 
also a discussion below) that H* should almost surely have absolutely continuous 
spectrum in 3D. Presumably, the correlation length truly diverges when one starts 
to approach the mobility edge and it will be extremely interesting to cover the 
missing case of E G [— E e (X), —E ]- 

This result can also be established by studying the density of states. Once the 
DOS is shown to be small below E* , the localization is a rather straightforward 
consequence of known methods (as say in [10] ). However, if one wants to study 
the behavior of the model at the closer vicinity of E* , this extra step can be an 
obstacle, as DOS increases. 



Unfortunately, in this work we are only able to descend to the values E* = 0(A 4 £ ), 
cannot claim that the correlation length indeed diverges. 
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As was pointed to us by the referee, the numerical results seem to qualitatively 
agree with the suggestion that the mobility edge is near —CX 2 (e.g. [17]). If one 
uses Cauchy variables instead of the box distribution, it is possible to calculate 
certain quantities explicitly for such variables: E(/(w)) = f(i) for a function / 
having a bounded analytic continuation to the upper half plane. In particular, 
DOS is computable, and it is not small for negative energies (namely 0(y/X)). 
The lower mobility edge in d = 3 appears to be positive on the basis of numerical 
studies (e.g. [TT], Figure 5), which suggests that the existence of high moments of 
the distribution of the potential plays a crucial role in the analysis. 

1.3. Major steps in the proof. In the unpublished notes, T. Spencer [18] pro- 
posed to prove the localization near the band edge using the multiscale analysis, 
with the initial volume estimates coming from the fact that the the density of states 
in the Lifshitz tail regime is small. To control the density of states he suggested 
to truncate the resolvent expansion at some optimal point. The corresponding 
Feynman graphs become superficially convergent after the suitable renormaliza- 
tion. In this paper, we complete the proof of the result announced in [18| . com- 
bining Spencer's perturbative approach with the Aizenman-Molchanov fractional 
moment method and developing the detailed estimates on the error terms in the 
renormalized expansion. 

The following representation for a Green function R(x, y) will be useful: 

Lemma 1.1. For any integer N and energies E that satisfy |j.^[ ) we have the 
decomposition 

AT-l 

R(x,y) = ^2A n {x,y) + ^ A N (x, z)R(z,y) , (1.10) 

n=0 26Z 3 

with Ao(x, y) = R r (x,y), and where the (real valued) kernels A n , An satisfy bounds 
E(A n (x,y)) 2 < {An)\E* (c(E*) -=\ e"V — l«-»l , n>l; (1.11) 

E\A N (x,y)\ < y^Ny. (C(E*) —=\ e-v-irk-i/l , jv">1; (1.12) 

where C{E*) = K\n E* for some generic constant K. 
The zero order contribution Aq satisfies 

< A ( X ,y)= f e^-y^ ff < j ^— - (1.13) 

J T 3 e(p)+E* (\x-y\ + l) 

for all x,y G Z 3 , and behaves asymptotically as 

Mx, y) = (l + 0(yW) + 0{\x- v n) x _ . (1.14) 



Lastly, we have 



MMx,y)f < KX * '' e -2V2W\ x - vl _ (L15) 

\x -y\ + l 

One then looks for the optimal value N to stop the expansion - note that the in- 
creasing factor of (4/V)! in An(x, y) competes with the decreasing factor (\ A E*) N / 2 . 
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The choice E* > A 4 ~ e has the effect that 

\2 

C{E*) -j= < X Bl , 0<B<1, (1.16) 

V E* 

which suffices to control (jl.lip - (|1.12|) . It turns out that the appropriate choice 
for N should satisfy 

(see the next section for details). In terms of the A - dependence, it corresponds 
to N ~ \~ be for b < B. Note that the square root in the denominator of (|1.16| is 
absolutely crucial for the strategy. To this end, let us mention that the represen- 
tation ()1.10p is a resolvent type expansion (cf. Lemma 13.11 below) . If one applies 
the rough norm bound on a each factor of the resolvent there, the denominator in 
(I1.16P will contain E* rather than its square root. The improvement is achieved 
using the Feynman diagramatic technique (Section [4])- 

Let us denote by H A,X the natural restriction of H* to A C Z 3 and let R A = 
(H AX + E + iO)^ 1 be the corresponding resolvent. 

Theorems Q] and [5] follow from the result above and from Aizenman-Molchanov 
a-priori bound on the fractional moment of the Green function as it appears in 
Lemma 2.1 of 3, which states that 

E\R A (x,y)\ s < C s (1.17) 

for any < s < 1, uniformly in x, y G A and A. Moreover, the bound above holds 
uniformly for an arbitrary set A. 

The rest of the paper is organized as follows: To make the presentation less 
obscure, we postpone the rather lengthy proof of the main technical Lemma 11.11 
until Section 2] and establish first Theorems [T] and [21 In Section [3] we perform a self 
energy renormalization required to get rid of the so called tadpole contributions. A 
technical statement regarding the properties of the self energy term a(E) is proven 
in Appendix [A"l 
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2. Proofs of Theorem Q] and [2] 



2.1. Theorem [TJ. Lemma ll.ll Aizenman-Molchanov a- priori bound ()1.17p . and 
Holder inequality imply that for any < s < 1/2 



E\R{x,y) - Rr(x,y)\ s 



N-l 



1 = 1 
N-l 



z£Z 3 



1/2 



|2 S \l/2 



< ^EAffoy) + ]T (Elljv^^l 28 ] ' (E| J R(z,y)| 2s ) J 

2 S \l/2 



i: ^(E^(x,y)) s/i + J2 (^\An(x,z)\) (E\R(z,y)\ 

z£Z 3 
K s/2 X s 



1 = 1 



< 



\x- y\ s / 2 + 1 
£((40!) s/2 (iT) 



JV-l 



/2B^|x-3/| 



■_, f C(E*)\ 2 y -.yffi\ X -y\ 



1=2 



<E* 



rv-i((4.Y)!r- a f^^V^ ^ e -V«i*-*i 



V \/IF / 

X S / 2 A S 
< e 

~ \x-y\ s / 2 + 1 

N-l 



z£Z 3 
-s\/2E' \x-y\ 



+ (C(E*)X 2 ) s e- s y^^^ J2 (( 4 00 



/2 /ccn^ 



i=2 



2\ 8(«-2)/2 



/J5* 



g(f) 



Choosing 



(47V) 4 



■((4J\0!) 



/£■* 



n s /2 



C(£*)A 2 \ siV/2 



/ 75* 



C(£*)A 2 



(2.1) 



one obtains, using the Stirling's approximation, that the summation over the index 
I is bounded by some s - dependent constant. On the other hand, for such N we 

have {my. ( C{ ^ 2 



N 



< e N . Hence, for such a value of N we have 



E\R(x,y) - Rr(x,y)\' < C s 



y e -sV2E* \x-y\ 

\x-y\ s / 2 + 1 



+ (C(B*)A a )*e-V¥l»-vl + (£*)-3/2 



exp 



C(£*)A 2 



(2.2) 



with some generic constant C s . We infer that for any e > and sufficiently small 
Ao(e) we have for any A < A (e) and E* > A 4 ~ e 



\\R(x,y) - Rr(x,y)\ s < C s 



A s 



\x- y\ s / 2 + l 



, \x-y\<{E*) 



-1/2 



(2.3) 
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(that is we proved Theorem [1]) and 

E\R(x,y)\ s < C' s (\ s e - s ^P\ x - y \ + e~ sX ~' /9 ) , (2.4) 

for any x,y £ Z 3 , where we have used the bound (|1.14| for R r (x,y). The latter 
estimate will be used in the proof of Theorem [2] 

2.2. Theorem [2l Using notation introduced after (J1.17J1 we define the decoupled 
Hamiltonian H\ to be: 

where A is a cubic box of the linear size 2L centered around the origin, and A c := 
I? \ A. Let dk be a boundary A. We will denote R A := {H A + E*)~ l . For any 
s < 1/2 and n £ dA we have 

E|i? A (n,0)| s < E\R{n,0)\ s + E\R A (n,0) - R{n,0)\ s 

= E\R(n,0)\ s + E\(R A (H A -H*)R)(n,0)\ s 

< E\R(n,0)\ s 

+ ]T {E|(i? A (ff A -^))(n,fc)| 2s } 1/2 {E| J R(fc ! 0)| 2s } 1/2 , (2.5) 

dist(fc,3A)<l 

where we used locality of (H\ — H£) - its non-zero matrix elements lies essentially 
on the boundary of A. Similar considerations lead to the estimate 

E\(R A (H A -H^))(n,k)\ 2s < C sup E \R A (n,m)\ 2s . (2.6) 

m: dist(m,<9A)<l 

Using bounds (|2.6p and (|1.17p as well as the Holder inequality, we obtain from (|2.5p 
E\R A (n,0)\ s < C S J2 {E|i?(fc,0)| 2s } . (2.7) 

fcGdA 

Plugging the bound (|2.4|) into the latter equation, with s < 1/4, we establish 

E|i? A (n,0)| s < C s L 2 (\ s e - sI VW + e~ sA " /9 ) (2.8) 

Now we are in a position to use the fractional moment criterion [3], Theorem 1.2, 
which states that if 

B s L 4 \- 2s ^ E\R A (n,0)\ s < b, (2.9) 

n(£d A 

for a certain constant B s and b < 1, then we have a bound 

E\R(x,y)\ s < lL- e W\*-v\. (2.10) 



It is easy to see from (|2.8[) that (|2.9[) is satisfied, provided that E* > A 4 c , and A 
is small enoughj, with 

V sVe* 

hence the result. 



Since the bottom of the spectrum of H^ is almost surely located at Xa, we can assume without 
loss of generality that E* < |a|A. 
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3. Renormalization OF TADPOLE'S CONTRIBUTION 

3.1. Expectation of the product of the random potentials. In what follows, 
the expectation of the product of random variable V u , namely 

n 
3=1 

will play an important role. These products naturally arise when one starts to 
expand the operator R := (H* + E + i0) _1 in the resolvent series 

n 

R = Y / (- XR oVu) l Ro + {-XRoV u ) n+l R 

t=0 

about the (unperturbed) operator R n :— (H + E)^ 1 . Indeed, 



(-XRoV u ) n Ro(xo,x n+1 ) = (-X) n Y. YlV^x^YlMx^Xi+x). (3.1) 

Xj£Z 3 ;j=l,..n j=l i=0 

Let T N , N < be the set {1, ..., N, N + 2, ...TV + N' + 1}, while IIjv.jv' will denote the 
set of partitions of Tn,n< into disjoint subsets Sj of size \Sj\ S 2N. Two partitions 
■k = {Sj}JL 1 , it' = {Sj}jL 1 are equivalent, n = n', if they coincide up to the 
permutation. For S C Tjf,N't let 

W=Ellt»h ( 3 - 2 ) 

yel, 3 j£S 

where 8 X , x G Z is Kronecker delta function, and xs denotes the collection of 
{xi , i G S}. One has an identity (see e.g. [5] Section 3.1 for details) 



n Wi) 

jeTjv,jv 



A? 



m=l 7 r={Sj}Jl 1 3=1 



(3.3) 



with coefficients c 2 i < (cl) 2l+1 proviso (A1-A2), and c 2 = E(V*(x)) = 1. The set 
Sj in the partitions 7r e H^at' can be of the special type: If 



S,- ={*,* + !} 



(3.4) 



we will refer to it as a tadpole, or a gate set. The non zero order contributions in the 
resolvent expansion, associated with the tadpole free terms, arc sufficiently small 
in the energy range (|1.4[) . as opposed to the contributions containing the gates. 

3.2. Self energy renormalization. The purpose of renormalization is then to 
include the tadpole contributions into the propagator itself. In our case this can 
be established by subtracting from the unperturbed operator — ^A the self-energy 
term a(E), described in Section [TJ We decompose H* as 

H r + V, H r :=--A-(x{E), V := XV U + a{E) . 



Hi 



The corresponding resolvent expansion for R defined in (|1.6|) takes the form 



R = Y J (- R rVyR r + {-RrV) 



n+1 



R. 



(3.5) 



i=0 
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with Rr as in (|1.5|) . Note that for any x E Z 3 

cr{E) = X 2 R r {x,x). (3.6) 

In place of (|3.1|) we get 

n n 

(-XRrV) n Rr(xo,Xn+l) = JZ T[(-*MXj) - °{E)) \{ R r (x U X i+1 ) 

Xj£l?;j=l,..llj = l i=0 

(3.7) 
If we open the brackets in (|3.7p , we obtain 

^2 Rr(xo,x 1 )6(x 1 )R r (x 1 ,X2)9(x 2 )...R r (x n _ 1 ,x n )6(x n )R r (x n ,x n+1 ) 

0,ay6Z 3 ;j=l,..n 

where 0(x) is either —XV UJ (x), or —<r(E) (whenever 9{x) = —a(E) we will refer to 
it as a bullet). Since a(E) — 0(X 2 ) for all permissible values of E, see Appendix 
[XI one can unambiguously define the order I (in powers of A) of the particular 
contribution 

Rr(xo,Xl)9(xi)R r (xi,X2)0(x2)-.-Rr(x n ^i,X n )9(x n )Rf,(x n ,X n+ l) , 

(with i?j being either R r or R) according to the following rule: Each factor of cr(E) 
counts as 2, while appearance of the random potential counts as 1, and we add 
up all the exponents to get the order of the term. For instance, the order of the 
expression 

R r (xo,xi)a(E)R r (xi 1 X2)XV ul (x 2 )Rr{x2,x 3 )a(E)R(x3,X4) 

is 5. 

To handle the renormalization of tadpole contributions properly, we decide at 
which value of n to halt the expansion in (|3.5|) individually for each contribution 
according to the following rule (to which we will refer as a stopping rule): If the 
order of the kernel 

R r (x o ,xi)0{xi)R r {xi,X2)O{x2)..-R r (x n -i,x n )9(x n )R{x n ,x n+1 ) (3.8) 

reaches or exceeds a value N to be determined later on, we stop expanding this 
particular term. To illustrate this procedure we write down the expansion obtained 
in a case of N = 2: 

R = R r - R r a{E)R - {XRrV u R} = 

R r — R r <r\E)R — XRj-V^jRj, 

+ XR r V UJ R r (j{E)R + X 2 R r V u R r V u R, (3.9) 

where the term in the curled brackets is the one we expanded according to the 
stopping rule. Note that the penultimate term is of order 3. It is not difficult to 
see that for a general N we get: 

Lemma 3.1. For any integer N we have a decomposition (used in Lemma \l.l\) 

N-l 

R(x,y) = Y^ (^2A , l (x,z)R r (z,y) + A , N {x,z)R(z,y) + B N (x,z)R(z,y)' S j 
zez 3 i=o 

N-l 

= J2Mx,y) + J2 A N (x,z)R(z,y), (3.10) 

1=0 z£Z 3 
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where A'^x, z) — 5\ x — z \, A'^x, z) is a summation over all possible terms of the type 
53 Rr{x,xi)0(x 1 )R r {xi,X2)9(x 2 )...R r (x n ,z)6(z) (3.11) 

6»,XjGZ 3 ; j=l,..n 

which are of the order I > ; while 

B N (x,z) = -a(E) 53 A' N _ 1 (x,w)R r (w,z). (3.12) 

wez 3 

The quantities Ai and An are defined as 

Ai{x,y)= 53 A' l (x,z)R r (z,y), A N (x,y) = A' N (x,y) + B N (x,y) . 

zeZ 3 

Proof. Note first that it is follows from (l3lT| l3~T2"l) that 

A N+X (x,y) = B N (x,y) - A 53 A' N (x,z)R r (z,w)V 0J (w)R(w,y) . (3.13) 

z,weZ 3 

We now prove p.lOp by induction: The base of induction, N = 0, 1, gives equalities 

R(x,y) = R(x,y), 

R(z,y) = Rr(z,y) - 53 (XR r (z,w)V LU (w)R(w,y)+aR r (z,w)R(w,y)). 
wez 3 

Suppose that (|3. 10|) holds for N, then 

w-i 
R(x,y) = ^2^^2A[{x,z)R r {z,y) + A' N {x 1 z)R(z,y) + B N ( : x 1 z)R{z,y) 
ze z 3 1=0 

JV-1 



X] ( X! A 'i( x > z ) R ri z i v) + A ' N (x, z)R r {z, y) + B N (x, z)R(z, y) 
ze z 3 1=0 

- 53 A' N (x,z)(\Rr(z,w)V u (w)R(w,y) +crR r (z,w)R(w,y) 



z,weZ 3 
N 



= 53 {^A[{x 1 z)R r {z,y) + A' N+l {x,z)R{z,y) + B N+1 {x,z)R{z,y)y (3.14) 
zez 3 i=o 
where in the last line we used (|3.12[ I3.13J) . □ 

Such a stopping procedure guarantees cancelation of tadpoles and bullets in the 
following sense: 

Lemma 3.2. We have 
i 
EA 2 (x,y) = X 21 53 53 53 -Rr(a;,a;i)i?r(a;,a;i+2) Jj R r (x l ,x l+ i) 
m=1 {-Sj}™ i (z 3 ) 2 ' »£f t ,, 

ro 
X n^l^K): ( 3 - 15 ) 

3=1 

wii/i S(z 3 ) 2 ' standing for summation over H? of all variables Xj with j £ T;.;, 
.T; + i = X2(+2 = y, and where ^2 denotes summation over all possible partitions of 
Tu which do not contain gates. 
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Equivalently, in the momentum representation 
EAf(x,y) = A 



2f„ ,.\ _. \2I I in d Pl+l dp 2 l+2 



e(pi +1 )+E* e( P2l+2 ) + E* 



dp t 



i 



>< n ^fr^E e n™*a>-*»o, (3.16) 

wftere 

a := 27r{-(pi +p J+2 ) ■ x + {p i+1 + P21+2) • y} ■ 

Remark: In order to obtain (|3.16|) from (|3.15p one uses 

J2tt(z-w)p " P __ / i2ir(z-w)p " P 



Rr(z,w) = / e * Ar i*- u '»- r rr = e l 

Proof. Let us introduce some extra notation. Let Ily;/. denote the set of partitions 
of Tj j into disjoint subsets Sj such that each of these partitions contains exactly 
k tadpoles (that is there are k subsets Sj of the form Sj = {i,i + 1}). For every 
partition LI; ^ . ^ 3 ir^ = {Sj}"L 1 we will denote by T(7Tfc) a subset of T^j which 
consists of the indices that differ from the set T c (7Tfe) := T; j ;\T(7r&) of the gate's 
indices associated with it^ (so card(T(7Tfc)) = 21 — 2k). Let itk = {Sj}sj.,s j G't(nk)}- 
Let v(y) be the connected segment of the set of the indices T c (7Tfe) which contains 
index v, that is 

v e i/(v) C T c (^ fc ) 

and 1/(1;) = {i,i + 1, ..., ft- — 1, ft}, with i — 1, ft + 1 e T(-7Tfc). Let d(v) denote the 
position of index v with respect to f(v), for example, if f (4) = {3, 4, 5, 6, 7, 8}, then 
d(4) = 2. We define iltieT ^ ) to be a product over such v € T c (7Tfc) that d[v) 
mod 2 = 1. For instance, if T c (7r/j) = {1,2, 4, 5, 6, 7}, the product will run over the 
variables 1,4,6. 

We can now express Ai(x,y) 2 as 

l 

Af( X ,y) = J2^~ 2k (-°( E )) k E E Rrfr*l)M*,*l+2) 

k=0 7r fc Gn,, i;fc ( Z 3)2l 

X II VufeJRriXi^i+i) Yl Rr{Xj + i,Xj +2 )5\ Xj + l - X] \. (3.17) 

ieT(7r fc ) jeT c (7r fc ) 

Note that here the index k corresponds to the number of bullets in the corresponding 
contribution, and is not related to the number of the tadpoles (which will show up 
as index k' below once we undertake the expectation over disorder). 
On the other hand, since C2 = 1, one obtains an identity 

N m 

E E Ih^K-) 

N N-k' m 

= EE E Lh^K) II W—I- ( 3 - 18 ) 
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At this point, let us first compute the expectation of k = part of the summation 
in (EU71) : 



iGT(Tro) 



A 2 ' Y^ Y Rr{x,Xl)Rr(x,Xl+2)E 

7roen iii; o (z3) 21 

I l—k m 

= X2L Y Y Y Y Rr(x,xi)R r (x,xi +2 )Y\c\ S] \8{x S] ) 

fe'=0m=l 7r k ,={S 3 }J =1 (23)2! j=l 

X Y\ Rr(Xi,X i+ l) Yl S \x v+1 -x v \ 

ier(i,i) verc(TT k ,) 

i i-k' 

= J2 X 2l ~ 2k ' (a(E)f J2 Y Y Rr(*,*i)M*>*i+*) 

fc'=0 m=l TTfc/K^}^! (Z3) a 

X II C l S 'jl <S ^' S 'j) II Rr{x l ,X l+1 ) Y[ Rr(Xj+l,X j+2 )5\ Xj+1 _ Xj \, (3.19) 

3=1 ier(7r k ,) iex-(7T fc ,) 

2 r? /„ _\ _ ^./77i\ r„^ „ii _ ,- nti 



where we have used X R r (z, z) = a(E) for all z e Z 
More generally, we have an equality 

A 2 '- 2fc (-a(£)) fc £ £ R r {x, Xl )R r {x,x l+2 ) 

7Tfc6n i-i;fc ( Z 3)2I 



x E 



Y\ V ul {x l )R r {x l ,x l+1 
ier(7T k ) 

l-k 



Yl Rr{x ] + l,X :j+ 2)5\ x . + 1 _ Xj \ 
jeT<=(7r fe ) 



(-cr(25)) fc (er(£))* 



„ , k'+k - 

= £(,,F 

fc'=0 

I-fe-fe' 

X Y Y Y Rr(x,X 1 )R r (x,Xi +2 ) 

m=1 *W+k={Sj}f =1 (I 3 ) 21 
m A 

X n !^! 5 ^' 5 '^ II Rr(xi,X l+1 ) Yl Rr(Xj+l,X J+2 )S lXj + 1 ^ X]l . 

i =1 ier(7r k , +k ) jerH^k'+k) 

(3.20) 

Therefore, the expectation of the rhs of (|3.17p is given by the formula 



EAf(x,y) = ]rA 2 <- 2 ^(£) £ ("^(T 

/3=0 fc,fc':{fc+fc'=,9} 

l-fi ^ __^ 

X Y Y Y Rr(x,Xi)R r {x,Xi +2 ) 

m=l ^ = {S 3 }™ 1 (Z3) 2! 
m A 

X II !^!' 5 ^ 5 '^ II Rr(x l ,X l+1 ) Y[ R r(%j+l,Xj+2)5\x j+1 -x j \, (3-21) 
3=1 ieT(^) j£T=(7r 3 ) 
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and the only non vanishing contribution comes from (3 = 0, since 

E ^) fe '(T) = (1 - 1} "' 

fc,fe':{fc+fe'=/3} 

hence ([3J5]l . D 

At this point we have to introduce additional notation (borrowed from [BJ): 

Definition 1. We consider products of delta functions with arguments that are 
linear combinations of the momenta {pi,p 2 , ■ ■ ■ ,P2n+2}- Two products of such 
delta functions are called equivalent if they determine the same afiine subspace of 
T 2n + 2 = {p t ,p 2 ,...,p 2n+2 }. 

One can obtain new delta functions from the given ones, by taking linear combi- 
nations of their arguments. In particular, we can obtain identifications of momenta. 

Definition 2. The product of delta functions A^ forces a new delta function 
6(^2 -j a jPj)i if ^2j a jPj — is an identity in the afiine subspace determined by 

A w . 

One can readily see that in the integrand of rhs of (|3. 16[) one has a forced delta 
function 5(pi - pi+i + pi+2 -P21+2), hence 

„ 2n+2 1 

A n>E .(x-y):=EA 2 n (x,y) = A 2 " J e -«*fei+p»+ a Mx-v) J| dpt _ _ 

n m 

x E E n%i^Eft-P'+i)- ( 3 - 22 ) 

4. Proof of lemma [TTT1 



Let us start with a remark that although the proof below is inspired (and closely 
follows) by work of Erdos and Yau [BJ, the geometry of energy surfaces plays no 
significant role here. In [IB], the geometry of energy surfaces poses a central difficult 
problem because the energy parameter varies, and in particular assumes values in 
the bulk of the essential spectrum of the nearest neighbor Laplacian. However, 
in the situation discussed in this paper, the reference energy E* is fixed, and lies 
below inf(cr(-A)) = 0. 

4.1. Reduction to the pairing case. In order to estimate A Ut E* (x — y) it suffices 
(up to the combinatorial factors) to consider the special case of partitions tt that 
appears in (|3.15p . where tt = {Sj} T j =1 with card Sj — 2 for each j - the so called 
pairing case. All other contributions are dominated by the corresponding pairing 
counterparts, as becomes transparent from the positivity of the free lattice Green 
function R r (x,y), cf. (|1.13p . Indeed, for any partition tt = {Sj}"L 1 choose an 
arbitrary subpartition tt' = {S'A" l =1 into pairs. Evidently, 

m 

y^ R r (x,xi)R r (x,xi +2 ) Yi Rr{xi,x i+ i) Y\_5{x S] ) 

(Z3) 2 ™ ieTi.i 3=1 

III: 

< ^ Rr(x,Xi)Rr(x,Xi +2 ) JJ Rr(Xi,X i+ i) J| 5(x S >.) , (4.1) 
(Z3) 2 " »6T l( i j=l 
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while the factor Y\7=i \ c \Sj\\ ^ s bounded by (cn) 2n+1 , see discussion in Subsection 
13.11 Hence, if one gets some bound M on the pairing type contributions, the whole 
A-n t E*{x — y) term can be rudely estimated as (2cn 2 ) 2n+1 M (where we took into 
the account the number of the possible partitions). 

4.2. Feynman graphs. A n: E-{x — y) is conveniently interpreted in terms of the 
so called Feynman graphs (the pseudograph, to be precise, since loops and multiple 
edges are allowed here). The graph, associated with particular partition 7r of T„ in 
is constructed according to the following rules (see Figure [1] and [2]) : We first draw 
two line segments, each containing n vertices (elements of T nj „). The vertices are 
joined by directed edges (momentum lines) representing momenta: p\ , . . . , p n +\ 
and Pn+2, ■ ■ ■ ,P2n+2- To each line pj we assign a propagator F(pj), with some 
given function F, save momentum lines p\ and p n +2, which carry additional phases 
e -i2„ Pl .(x-y) and e -i27rp„ +2 -0- y ) 5 respectively. For tt = {Sj}f =1 we identify all 
vertices in each subset Sj as the same vertex (in Figure [TJ the paired vertices are 
connected by dashed lines). Note that thanks to the existence of the forced delta 



p 4 



Pi 



Pio 



Figure 1. Construction of the Feynman graph, part I, n = 4. The 
corresponding delta functions are §(p\— P2+P3— Pa), &{Pi—P5+P9 — 
Pio), 5(jp2—P3+P&— Pt), and S(p7~pg). The last delta corresponds 
to the tadpole. Note that the sum of all momenta in the above delta 
functions gives a forced delta function S(pi —P5+P6 —pio), hence 
we can introduce the dashed lines connecting vertices 1,6,7, and 
12, identifying them as a single vertex. 



function 5{p\ — pi+\ + P1+2 — P21+2), we can identify vertices 1, n, n + 1, 2n as a 
single one, and therefore one can think about the closed graph (with special rules 
that apply for momentum lines p\ and p n +2, mentioned above). To summarize, 
the outcome of this construction is a directed closed graph, which is called the 
Feynman graph associated with the partition tt. The momenta in the graph satisfy 
the Kirchhoff 's first law, that is the total momenta entering into each internal vertex 
add up to zero (if arrow faces outward the vertex, we count its momentum with a 
minus sign). A tadpole corresponds to the so-called 0-loop, that is some (directed) 
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line of the graph claims one vertex as its both endpoints. For a given Feynman 
graph G, one can choose a particularly useful expression for the product of delta 
functions A„. Choose any spanning tree of G which does not contain momentum 
lines pi,p n +2- The edges belonging to the spanning tree will be called the tree 
edges (momentum lines), and all the rest are the loop edges (since an addendum of 
any loop's momentum line creates a loop). Let us enumerate the tree variables as 
Mi,...,Ufe, and loop variables as Wi,...,wi, with say w\ =pi,ui2 — Pn+2 (note that 
k + I = 2n + 2). The number k of the tree momenta coincides with the number of 
the delta functions in A w . One can check (see e.g. [B]) that the product of delta 




Figure 2. Construction of the Feynman graph, part II: Identifi- 
cation of the vertices. The tadpole corresponds here to 0-loop. 

functions A w is equivalent to 

k I 

Yldim-Y^aijWj), (4.2) 

i=\ j=l 

with 

{±1 loop is created by adding Wj to the spanning tree contains Ui 
otherwise 

The choice of the sign depends on the mutual orientation of Ui and Wj . 

4.3. Exponential decay. We first want to establish the exponential decay of 
A n ,E*{x) that appears (|3.22[) from the simple analytic argument, and then in the 
next section obtain the bound on A„^;»(0). 

Specifically, we want to show that for a general value of n, 



A n ^(x)<e-^^^A rhE , /2 (0). (4.3) 

Indeed, for any given x £ 1? lets choose 7 G {1, 2, 3} such that 

\x ■ e-»| — max Ix-eJ. (4-4) 

ie{l,2,3} 
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Then \x ■ e 7 | > |x|/\/3- In order to obtain the exponential decay of A n: E*{x) we 
first perform the integration in the rhs of (|3.22[) over the tree momenta, using 
Let us use the shorthand notation ^^ for a sum over all possible partitions in 
c w for a product of the corresponding csj, and r v will denote the number of the 
delta functions containing the loop momentum w\ in the tt's partition. We get 



— i2-KW\-x 



2n + 2 



Ju"p"- a — "if 7^. <«) 

J te*' i=rv + l W 

where $' is a set of all loop variables in the partition ir, except for wi. The variable 
qk is some linear combination of the loop variables in $', for k = 1, ..., 2n + 2. We 
want to obtain a bound on the integral over wi variable. Note that 



*n 



e 



-i27rp-:r 



L = i e(p + qt) + E* 

1/2 



= f dp' e- i27ripx - pe -' xe -' ) [ ' d{p ■ e 7 )TT -7- ^ — e -^( P - e ^-e T ) 

J J -111 ;_, 



(4.6) 



where J dp' stands for integration over components of p orthogonal to e 7 . Without 
loss of generality, let us assume that x-e-y > 0. It is easy to check that the integrand 
as a function of p-e 7 is 1-periodic, analytic inside the rectangle formed by the points 

{-1/2; -l/2 + n/lF/5; 1/2 + n/iF/5; 1/2} 
for sufficiently small E* : Indeed, we have 

Re e(p + q + iee 7 ) + E* > e(p + q) + E* /2 
uniformly in q, provided 

Ve* 

0<e< — — 



for small E* , where we have used sin(a + ib) = sinacosh6 + ?cosasinh6. Moreover, 
the periodicity implies that the integrals over the vertical segments coincide: 

-1/2+i*/E v /(2tt) J__ -, 

— i27r(p-e- f x-e-„) 



/-J./2+JViS*/(.27TJ 
d(p ■ e 7 ) TT 
-1/2 ti 



1-1/2 1J ^\e(p + qi )+E* 

fl/2+iVW/(2n) r 
/ dip • e 7 ) TT e -i27r(p.e 7 x.e 7 ) u .7) 

A/2 l 7; fJ e(p + %) + £* 
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Therefore 

,1/2 



dip ■ e 7 ) TT "T - > i f* 

1/2 T=\ e (P+<l')+ E 



j—i27r(p-e- / x-e- / ) 



l/2+i\/W/(2ir) 



d (p • e 7) n 



1/2 

•: e 



e(p + ft) + E* 



-i27r(p-e-yX-e-y) 



^^ I' d(p-e y )f[ 

J -V2 i=l 



1 



1 ' _I_ l 



< e -MVB=73 r d(p-e 7 )n 



eb + *)+^*/2 



(4.8) 



hence 



4.4. Bound on the pairing type diagrams. Note that for any p 6 T 3 we have 
an elementary estimate 

e(p) = 2^sin 2 (^. ei ) > 2 sin 2 f*! 



(=i 



Vy^ 



On the other hand, by Jordan's inequality sin 2 I ^& ) > -|- for any p £ T 3 . There- 
fore, 

{e{p)+E*)- 1 < C{p 2 +E*)- 1 (4.9) 

for any p£T 3 , and 



2n+2 



n * 



c 



pf + E* 



m=l *={•%}£,! \i=i / ^Si 

Now consider the pairing type partition n. For such a partition all c\g.\ = 1. It is 
natural to rescale variables, p t = y/Wqt , to get 



2n+2 



c 



n ^« 2 , E * sc^Pi-Pi+i) 






/' 2n+2 „ n 



where we used the scaling S(ap) = a~ 3 ^(p) in three dimensions (and the fact that 
there are n delta functions involved in the pairing case) . Let us note that n is equal 
to the number of loop momenta and equal to the number of vertices minus 1. Each 
integration now runs over the torus (i?*)~ 1 / 2 T 3 , concisely denoted as J . We can 
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only increase the right hand side by replacing it with 



ieSi 



l n 8n+8 £* r 2n+2 



(^ ) 2 V 3 ) 2n+2 f=i fe + 1 ) 111 (% + 2 ) J= i ,t^ 

(4.12) 

where <$(•) in the last line stands for the standard Dirac delta distribution on K 3 . 
At this point we are going to switch to the (Euclidean) Feynman graph repre- 
sentation, as described in Subsection 14.21 with the propagator 

F{qj):= (g J 2 + l)ln 4 (,| + 2)- 
Let us define the value \G\ of the graph G (characterized by its partition -k = 

/2n+2 1 n 

R3)2 „ + 2 ax (g£ + 1) In (# + 2) ax z_. 

If partition 7r is tadpole-free (and those are only partitions that enter into (|3.22p 
thanks to the self energy renormalization) , the corresponding Feynman diagram 
has an important property of being (logarithmically) superficially convergent. The 
value of such graphs is controlled by the following theorem (an adaptation of The- 
orem A\ in [7] to the case in hand): 

Theorem 3 (Bound on superficially convergent Euclidean Feynman graphs, [7]). 
If Feynman graph satisfies assumptions A\ -AS below, then its value is bounded by 
K n , with some generic constant K . 

This result yields that the contribution from a pairing type partition to A n ^E* (0) 

is bounded by 

/ A 2 \" 

E* [C \n 9 E* -£=) (4.14) 

- the key estimate of this subsection. 

The subsequent assumptions require the introduction of some additional notation. 
Let a subgraph G' C G be a subset of the lines of G. The vertices of G' are the end 
points of lines of G' and an external line of G' is an edge of G \ G' which is hooked 
to a vertex of G' . 

Let us denote by A(G') the number of loop edges of G' (defined in Subsection 
I4.2p . by E(G') the number of external edges of G' , by I(G') the number of internal 
lines of G' , and by N(G') the number of vertices in G' . 

Two subgraphs F, F' are disjoint if they have no line and no vertex in common. 
They overlap if they are not disjoint and do not satisfy an inclusion relation (F C F' 
or F' C F). A forest F is a set of non-overlapping connected subgraphs. 

A subgraphs F is said to be one line reducible (OLR) if there exists a line of 
F such that its removal increases the number of connected components of F. A 
connected subgraph which is not OLR is called proper. 
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.41 The lines 

Each line has a propagator of the form 

1 



( 9 2 + l)lnV + 2)' 
A2 Superficial convergence 

For any connected subgraph G' C G a superficial degree of divergence 
div(G') is given by 

div(G) = 3A(G) - 27(G) 

For any connected subgraph G' C G a logarithmical superficial degree 
of divergence / — div(G') is given by 

l-div(G') = A(G') - 4/(G'). 

A graph G is called superficially convergent if any connected subgraph G' 
of G satisfies either 

div{G') < -2eE{G') 

or 

dw(G') = and I - div(G') < -e . 
Ai div — Forests 

(i) There exists a constant G such that the number of proper div = forests 
of G (i.e. forests consisting of proper subgraphs G' with div{G') = 0) is 
bounded by G L(G) . 

(ii) There exists a constant G' such that for every connected G' C G with 
div(G') — 0, G' has at most G' external vertices. 
Remark: Theorem 41 in jJJ is a much more general result than the one pre- 
sented here. It contains two additional assumptions which are irrelevant in our 
context (namely HA.2 — HA.3 in [7]). We also adapted the various definitions from 
[7] (such as a superficial degree of divergence) to the concrete situation discussed 
in this paper. 

The rest of this subsection is devoted to the validation of the assumptions A2-A3 
of Theorem [3] in the present context. A3 is met since all superficially divergent sub- 
graphs in our situation (in fact there is only one divergent subgraph F, introduced 
below on Fig. [3]) are also divergent in <f>\ theory, where this assumption holds true 
with G = 8 (4J). The constant C' (which in our context corresponds to the number 
of external lines of the aforementioned graph F) is equal to 2. 
We now want to establish the validity of A2, with e = j^: 

Let us note that for a pairing partition, the degree of each internal vertex is 
4, that is we have 4-regular directed graph. For such a graph any spanning tree 
contains n edges, and n + 2 loops accordingly. 

Since G is 4-regular, it is easy to see that for any connected subgraph G' 

I(G')< m ( G ')^ E ( G ') . A(G') + N(G')-1 = I(G'), (4.15) 

where the latter relation follows from the fact that the spanning tree for G' contains 
N(G') — 1 lines, and the rest of the internal lines can be thought of as loops. Hence 

div(G') = -A(G') + 2(2A(G') - I(G')) 

= -A(G')+2(I(G') + 2-2N(G')) < 4 - E(G') - A(G') (4.16) 
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and 
l-div{G') = (A(G')-I(G'))-3I(G') = 1-N(G'))-3I(G') < -4 (4.17) 

for any subgraph G' of G with N(G') > 2. Since N(G') = 1 corresponds to the 
tadpole, which is not allowed, we conclude than is that all relevant subgraphs are 
logarithmically convergent. 

For the whole graph G we have div(G) = 3(n + 2) - 2(n + 2) < -n/4 for n > 2, 
while div(G) = for n = 2. Since 2E(G') mod 4 = for any subgraph G' of the 
4-regular graph G, and #(G') ^ unless G" = G, we deduce from (|4.15p and (|4.16|) 
that for any Feynman graph G corresponding to (|4.12|) with n > 2 the only possible 
proper subgraphs G' with div(G') > are: 

(1) N(G') = 1, A(G') = 1, E(G') = 2 - a 0-loop, that is a tadpole. For the 
0-loop div(G') = 1, but on the other hand, the tadpoles are prohibited in 
our partition. 

(2) N(G') = 2, A(G') = 2, £'(G / ) = 2 - either a pair of tadpoles connected by 
an edge or a graph F, shown on Figure [3] (with two external edges omitted). 
For the latter graph we have div(F) = 0, hence F is superficially convergent 
as well. 




Figure 3. Graph F. 

Since it is clear from (|4~16ll that div(G') < -E(G')/5 for E(G') > 5 while div(G') < 
— 1 < —E(G')/h for E(G') < 5 (with the exception of the graph F), we conclude 
that all tadpole- free pairings are indeed superficially convergent. 



4.5. Bounds (|1.11|4TTT5|) . We deduce from (|4.14p . (|4.ip and discussion thereafter 
that 

A n>E *(0) < (in)\E* (c \n 9 E* -A 



and (|l.lip follows now from (|4.3p . 

To get (j!.12|) note that it follows from Lemma 13.11 that 

A N = A N (-^A + E*) - a(E)A N _ 1 R r . 
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We therefore obtain 



E\A N (x,y)\ < J2 {{E\A N (x,z)\ 2 ) 1/2 • \(-1a + E*)(z,y)\ 

z£Z 3 *• 

+ a(E) (E\A N ^( Xl z)\ 2 ) 1/2 ■ R r (z,y)} 
< 7 ]T (®\An(x,z)\ 2 ) 1/2 

z£Z 3 : \z-y\<1 

zez 3 ' ' 

with some generic constant C, provided E* < 1 and where we used the bound 
(|1.14p on the free Green function R r (z, y). 

It is clear from (jl.lip that the first contribution is bounded by 



W\ x ~y\ 



f 


A 2 ,N/2 

-/= e 

Ve* J 

-V2W\z-y\ 


C'^J{AN)\E* (C ln y £* 
lile the second one is bounded by 
n' p(\ f*\ V^ r — \p^\x-z\ 


< C'PiX^E^e-^-^-y 1 Y 

zez 3 


1 
1 i 7 e 

\z-y\ + \ 



-V2S*(l-l/V24)|z-j/| 



< ^(A.FJe-v^M 



with 



P{\, E*) := A 2 ^J{AN -4)\E* [C In 9 E 



A 2 N AT/2-1 



IE* 

hence (fl~T2|) . 

The positivity of Aq(x, y) in (|1 . 13|) is immediate from the positivity of e (cc, y) 
for all non negative values of i, which in turn is obvious, since all off diagonal matrix 
elements of the Laplacian are positive, and its diagonal part is proportional to the 
unity operator. 

The upper bound in (|1.13| is the standard estimate, see e.g. [12], it also follows 
from (11.14J) . The bound (|1.14|) for the free Green's function is known as well, see 
e.g. 0. 

The convenient way to establish bound (|1.15|) is in the coordinate representation: 
Since E(V LJ (x)V LJ (y)) — 5\ x - y \ one checks that 



EA\{x,y) = \ 2 ^i? 2 (z,z)i? 2 (z,y). 
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Now we use (|1.14|) to bound the rhs from above as 



A 2 V l 1 

^ \x- z\ 2 + 1 \y - z\ 2 + 1 



2V2E^(\x-z\ + \y-z\) 

y 



zei 

1 1 



< A 2 e -2V2^|x- B | ^ 



, x-z\ 2 + l \y-z\ 2 + l 

~ \x-y\ + l V ' 

Appendix A. Properties of the solution of (|l.3|) 
It is instructive to rewrite ljl.3|) as 

,2 / d 3 p 



E = E* +\ z / , s " . (A.l) 

J T 3 e(p) + £* 



S then is a well defined function of E* on the ray E* £ [0, oo), with 

E(0) = A 2 



X i / d3 P 



/is e(p) ' 
Differentiation of (jA.ip with respect to 25* gives 

E'(E*) = l-\ 2 [ , -, f*,. 2 ■ (A.2) 

J T 3 (e(p)+.E*) 2 

Monotonicity of the integral on the rhs of the above equation for E* E (0, oo) 
implies that there exists exactly one extremum of E (namely minimum) for such 
values of E* . In particular, the function E = E(E*) is invertible for E > E(0). 
For any E* > we have an estimate 

d3 P ^ „tT**s-X/2 



L (e(p)+E*r ~ C{E * 



which follows from (|4.9[) and the extension of the domain of integration to the whole 
E. 3 . We infer that 

E(E*)= f E'(t)dt + E(0)>\ 2 f ^P-+E*-C\ 2 f t~ l/2 dt, 
Jo Jjz e(p) J 

hence E(E*) > E(0) for E* > 4C 2 A 4 , so there exists an inverse function E* — 
E*{E) in this range. On the other hand, note that e(j>) < p 2 /2, hence 

d 3 p f d 3 P f d 3 p _ r , (F , ) -l/2 



T 3 (e(p) + E*) 2 ~ 7 T3 (p 2 + E*) 2 ~ J Be , (p 2 + E*) 2 

where Be* denotes the ball of radius VE* around the origin (E* is assumed to be 
small). Therefore we obtain 

„ f d 3 p <- E ' 



I - e(p) 
so that if 



<A 2 / 44+£* -C*'^ 2 / t'^dt, 
Jjs e(p) 7 



S >A 2 / ^ + A 4 -. 



/T3 e (p) 

then i?* > C*A 4 ~ e for small enough values of A, as follows from the solution of the 
corresponding quadratic equation for \/E*. 
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